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I.  INTRODUCTION 


During  the  last  two  decades  multi-body  dynamics  has  been  an 
active  research  subject  in  the  field  of  spacecraft  attitude  dynamics 
and  control,  and  in  the  field  of  mechanisms  and  manipulators.  Most 
of  this  research  has  been  oriented  toward  developing  computationally 
efficient  techniques  for  solving  (or  "simulating")  a set  of  nonlinear 
dynamics  equations.  Thus,  many  different  forms  of  the  equations  of 
motion  have  been  obtained,  and  studies  have  been  performed  on  the 
equivalences  and  relative  merits  of  various  alternative  formulations 
Cl]  - [3]  . All  of  the  various  formulations  can  be  shown  to  be 
equivalent  in  the  sense  that  they  can  be  obtained  from  each  other 
via  nonsingular  linear  transformations. 

The  purpose  of  this  paper  is  to  identify  some  aspects  of  linear 
and  nonlinear  system  theory  which  have  direct  application  to  multi- 
body dynamics.  Classical  dynamics  and  system  theory  have  common 
roots  in  the  transformation  theory  of  ordinary  differential  equations; 
consequently  it  is  not  surprising  that  many  system  theory  concepts 
have  applications  in  classical  dynamics  in  general,  and  in  multi- 
body dynamics  in  particular. 

This  paper  is  aimed  at  system  theorists  interested  in  applications 
to  classical  dynamics,  and  at  dynamicists  interested  in  system  theory 
techniques.  In  order  to  have  a common  background  for  discussion,  the 
paper  first  describes  the  relevant  aspects  of  a typical  spacecraft 
attitude  dynamics  and  control  system.  This  is  then  followed  by  a 
discussion  of  some  aspects  of  system  theory  which  can  meaningfully 
be  applied  to  the  spacecraft  attitude  dynamics  problem.  The  discussion 
is  in  part  informal  since  the  intent  is  to  demonstrate  applicability  of 
the  concepts  rather  than  to  prove  concrete  results.  The  subsequent 
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section  is  a description  of  multi  -body  dynamics  from  a vector  space 
point  of  view.  Constraints  on  a system  are  related  to  uncontrollability. 

Numerous  references  to  applicable  system  theory  are  provided 
for  the  interested  dynamicist,  and  some  key  multi-body  dynamics 
references  are  provided  for  the  interested  system  theorist. 

It  is  hoped  that  this  paper  will  contribute  to  the  cross- 
fertilization of  the  fields  of  system  theory  and  spacecraft  attitude 
dynamics.  Much  future  work  remains  to  be  done  in  this  area.  The 
problems  are  expected  to  continue  into  the  next  century  as  larger  and 
iarger  spacecraft  (some  several  miles  in  size)  are  being  considered 
and  studied. 

The  motivation  for  undertaking  the  study  reported  herein  came 
from  the  desire  to  find  computationally  efficient  spacecraft  attitude 
dynamics  formulations.  This  required  a systematic  study  of  the 
structure  of  various  alternative  formulations  [3].  As  demonstrated 
in  this  paper,  many  concepts  of  system  theory  are  natural  tools  in 
this  study. 


II.  SPACECRAFT  ATTITUDE  DYNAMICS  AND  CONTROL  SYSTEMS 


For  the  present  purposes,  a typical  spacecraft  attitude  dynamics 
and  control  system  can  be  considered  to  consist  of  three  parts: 
dynamics,  kinematics,  and  controller.  The  dynamics  and  kinematics 
together  form  the  plant;  this  plant  is  a "given"  which  the  control 
systems  engineer  normally  does  not  change.  The  task  then  is  to 
design  a feedback  controller  such  that  the  combined  system  of  plant 
and  controller  has  acceptable  performance.  A controller  consists 
of  sensors,  data  processors,  and  actuators;  these  components  perform 
the  functions  of  measurement,  computation  and  decision,  and  exerting 
force  (or  torque)  on  the  plant,  respectively. 

A.  Controller 

Controllers  are  always  nonlinear  because  real  physical  sensors 
and  actuators  are  always  nonlinear;  in  addition,  usually  it  is  also 
desireable  to  use  nonlinear,  and  perhaps  variable-structure,  data 
processing.  With  only  a small  risk  of  oversimplifying  matters,  it 
can  be  said  that  there  are  as  many  different  spacecraft  controllers  as 
there  are  different  spacecraft,  or  as  there  are  spacecraft  control 
systems  designers.  However,  often  a controller  can  be  modelled 
with  a nonlinear  finite -dimensional  state  equation,  e = 9 (e,  q,  t), 
and  a nonlinear  output  equation,  K = t(e,  q,  t);  here  e is  the 
controller  state  variable,  q is  the  plant  output,  K is  the  controller 
output  (force  or  torque),  and  t is  the  time.  Sometimes  a portion 
of  the  controller  is  discrete  rather  than  continuous.  Occasionally 
the  continuous  portion  of  the  controller  is  neglected  and  the  entire 
controller  is  modelled  as  discrete. 


B.  Dynamics 
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Whereas  the  form  of  the  controller  equations  are  at  the 
control  system  designer's  discretion,  the  form  of  the  plant  equations 
is  not.  However,  the  interesting  fact  is  that  the  plant  equations 
always  have  the  same  form  (except  at  various  discrete  times  the 
order  (dimension)  of  these  equations  may  change).  The  plant 
equations  can  be  written  as 

G = yG  + K (2.1a) 

g = VG  (2.1b) 

where  the  dynamics  state  variable,  G,  is  the  system  momentum, 
and  the  dynamics  output  variable,  g,  is  the  system  velocity.  G 

and  g are  elements  of  linear  vector  spaces,  and  therefore  they  are  ) 

linearly  related: 

G = M-g  (2.2) 

where  is  the  system  mass  (or  mass /inertia).  p.  is  always 
positive  definite  symmetric,  and  therefore  it  has  a positive 
definite  symmetric  inverse,  v = p.  . Since  g is  related  to  G via 
a nonsingular  linear  transformation,  it  follows  that  g itself  can  be 
used  as  the  dynamics  state  variable.  The  plant  equations  can  then 
be  written  as 

g = - yt  g + vK  (2.  3a) 

g = g (2.3b) 

where  y 1 is  the  transpose  of  y.  y is  called  the  system  Cartan 

operator  [4]  . Apparently,  Equations  (2.1)  and  (2.3)  are  formally  | 
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duals  or  adjoints  of  each  other  [5  ] - [6  ] , since  v is  symmetric 
( v = vt  ).  However,  the  definition  of  dual  or  adjoint  equations 
given  in  linear  system  theory  does  not  apply  since  equations  (2.1) 
and  (2.  3)  are  not  linear  equations.  The  nonlinearities  in  these 
equations  are  due  to  the  facts  that  y depends  (linearly)  on  G or  g, 
and  both  y and  v depend  (nonlinearly)  on  the  position  p . But 
equations  (2. 1)  and  (2.  3)  are  truly  duals  of  each  other  for  reasons 
whose  roots  lie  in  linear  algebra  (or  tensor  analysis):  G and  g 
represent  the  same  physical  quantity  in  different  vector  spaces 
which  are  the  duals  of  each  other.  Let  qy  he  the  velocity  space  so 
that  g e qy  ; then  the  dual  space  q/  is  the  momentum  space  so  that 

:,'c 

G f ^ . The  relationships  between  and  V can  also  be  expressed 

as 


qr  = 

p.qy 

(2.4a) 

qy  = 

vqy* 

(2.4b) 

Thus,  p maps  q/  into  its  dual,  and 

'*>  maps  q/# 

into  its  dual.  The 

vector  spaces  q/  and  q/  are  turned  into  inner  product  spaces  by 
defining  inner  products  as  follows  [7] 

< g.  g'>^  = gt  M-g'  (2.  5a) 

<G,  G'^*  = GtvC'  (2.5b) 

where  g , g'  ety  and  G,  G'  € qy  • But  the  kinetic  energy  T can 
be  expressed  as  [3]  - [4] 

T = \ G*g  = i-gVg  = j Ct  vG  (2.6) 

Combining  equations  (2.5)  and  (2.6)  now  yields 
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(2.7) 


2T  = <§’  g>qy  = <G’  G\* 

In  the  language  of  tensor  analysis  [8]  , the  momentum  G is  also 
referred  to  as  the  covariant  velocity,  and  the  velocity  g is  also 
referred  to  as  the  contravariant  momentum,  g.  is  the  covariant 
metric  tensor,  and  its  inverse  v is  the  contravariant  metric  tensor. 

The  force  K can  also  be  referred  to  as  the  covariant  force;  the 
contravariant  force  k is  defined  by 

k = v K or  K = \i  k (2.8) 

Evidently,  these  two  equations  are  analogous  to  equations  (2.1b)  and 
(2.  2).  The  time  derivative  of  the  kinetic  energy  can  now  be  expressed 
as  [3  ] 

T = Kg  = k*G  (2.  9) 

Note  that  T = 0,  or  T is  constant,  when  K = 0.  In  the  language  of 
tensor  analysis,  T and  T are  invariants  since  they  are  equal  to  the 
scalar  product  of  a contravariant  and  covariant  vector. 

The  Cartan  operator  y not  only  appears  in  the  time 
derivative  of  G and  g,  but  it  also  appears  in  the  time  derivatives 
of  p.  and  v as  follows  : 

m.  = yn  + ny1 

v = _ytv  - v y 

■ 

These  equations  are  useful  in  relating  equations  (2.1a)  and  (2.3a).  For 
; example,  differentiating  G = p.  g with  respect  to  time  yields 

■' 

! 


(2.10a) 

(2.10b) 
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G = P g + ng  I substituting  for  p from  (2.10a)  and  for  g from 
(2.3a)  then  yields  G as  in  (2.1a).  Equations  (2.10)  are  formally 
similar  to  the  linear  matrix  equations  which  appear  in  linear  system 
theory  (see  [5]  , pp.  58-62);  thus,  equation  (2.10b)  can  be  considered 
to  be  the  adjoint  equation  associated  with  (2.10a). 

All  the  dynamics  equations  given  thus  far  have  been  stated 
without  proof.  These  equations  can  be  derived  from  Newton's  law 
for  a particle  in  a system  of  particles.  However,  all  the  dynamics 

equations  can  easily  be  interpreted  by  using  the  tensor  analysis 

, , □ o 

concept  of  a covariant  derivative.  Let  G and  g be  defined  by 


□ 

G 


G - y G 


□ t 

g = g + y g 


(2.11a) 

(2.11b) 


and  let  p and  v be  defined  by 


a • t 

P = P - yp-py 


v=  v + y v + vy 


(2.12a) 
(2. 12b) 


A quantity  with  an  open  square  dot  on  top  will  be  called  the  covariant 
time  derivative  of  that  quantity.  Equations  (2.11)  and  (2.12)  follow 
from  the  covariant  differentiation  equations  of  tensor  analysis  when 
it  is  remembered  that  G is  covariant  of  degree  one,  g is  contra- 
variant  of  degree  one,  p is  covariant  of  degree  two,  and  v is 
contravariant  of  degree  two.  The  elements  of  y can  be  expanded 
linearly  in  terms  of  the  components  of  g,  the  coefficients  in  these 
expansions  being  the  components  of  the  affine  connection.  Using  these 
covariant  time  derivatives,  equations  (2.1a),  (2.3a),  and  (2.10) 
become  simply 
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and 


On 

ii 

K 

(2.13a) 

g = k 

(2.13b) 

a 

p = 0 

(2. 14a) 

□ 

v = 0 

(2.14b) 

Equations  (2.13)  state  simply  that  the  covariant  time  derivative  of 
the  covariant  velocity  is  equal  to  the  covariant  force,  and  the  covari- 
ant time  derivative  of  the  contravariant  velocity  is  equal  to  the  contra- 
variant  force.  Equations  (2. 14)  state  that  the  covariant  time  derivative 
of  the  metric  tensor  and  its  inverse  are  zero.  These  results  are  well 
known  in  tensor  analysis  [9  ] - [10  ] . Note  that  covariant  differentiation 

is  a derivation  [ll]  in  the  sense  that  it  follows  the  rule  of  Leibnitz: 

o o □ 

G = pg  implies  that  G = p g + pg  ; from  this  and  equations  (2.13) 
and  (2.14)  there  now  follows  that  K = pk  , in  agreement  with  (2.8). 

It  is  evident  from  Equations  (2.1)  to  (2.14)  that  the  dynamics 
equations  can  be  put  into  a very  neat  and  simple  form.  It  is,  of 
course,  also  possible  to  put  these  equations  into  more  complicated 
forms,  such  as  introducing  coordinates  with  respect  to  a time- 
varying  reference:  then  additional  terms  arise  in  the  expressions 
for  the  covariant  derivatives  and  in  the  equations  of  motion,  and 
the  kinetic  energy  is  no  longer  a quadratic  in  velocity  or  momentum. 
However,  the  resulting  more  complicated  equations  are  not  any 
more  general,  even  though  they  may  be  more  convenient  for  some 
particular  purpose. 

The  dynamics  equations  which  have  been  presented  are  not  in 
a form  which  are  familiar  to  most  dynamicists.  In  fact,  it  is 
usually  not  convenient  to  obtain  the  Cartan  operator  y as  a prelude 


-12- 


to  obtaining  the  equations  of  motion.  However,  if  y is  found,  it 

usually  can  be  decomposed  as  y = Vq  + y£-  where  y^,G  = 0 so 

that  y G = y G.  Similarly,  y can  be  decomposed  as  y = y + y ' 
b § § 

where  y'*  g = 0 so  that  ytg  = ylg.  Thus,  only  y„G  appears 
S & 

in  the  momentum  formulation  equation  of  motion,  and  only  y g 
appears  in  the  velocity  formulation  equation;  consequently,  it  is 
often  not  evident  that  the  same  Cartan  operator  y appears  in 
G and  g,  and  in  m>  and  v . 

A method  of  finding  y for  any  particular  problem  is  described 
in  Section  IV.  B.  Briefly,  the  method  consists  of  first  expressing 
the  equations  of  motion  in  terms  of  variables  for  which  y is  known  , 
and  then  transforming  to  more  convenient  variables.  Thus,  for 
example,  the  y for  a single  particle  (which  obeys  Newton's 
laws  of  motion)  is  zero,  and  the  y for  a system  of  particles  is 
therefore  also  zero.  By  considering  a rigid  (or  perhaps  flexible) 
body  to  consist  of  N particles,  the  y for  the  rigid  (or  flexible) 
body  can  be  found.  From  this  y the  y for  any  arbitrary  system 
of  bodies  can  be  found. 


C.  Kinematics 


The  system  momentum,  G,  and  system  velocity,  g,  are  vectors, 
but  the  system  position,  p , is  a point.  Thus,  G and  g are  elements 
of  a vector  space,  whereas  p is  an  element  of  a manifold  [ 12  ] - [ 13  H . 
If  the  plant  has  n degrees  of  freedom,  then  the  linear  vector  spaces 
A/  and  T;  have  dimension  n,  and  so  does  the  manifold  Therefore, 

g and  G have  n independent  components,  and  the  point  p is  specified 
via  n coordinates,  q^,  for  i = 1 to  n.  Let  q be  the  n-element  column 
matrix  which  has  q.  as  the  ith  element,  and  let  g now  represent  an 
n-element  column  matrix  whose  ith  element  is  the  ith  component  of 
the  system  velocity  vector  (with  respect  to  some  basis).  The  time 
derivative  of  q can  now  be  linearly  related  to  g: 
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q = ®g 


(2.15) 


where  a depends  (nonlinearly)  on  the  point:  O'  = a (p),  or 
a = a (q).  Note  that  q is  an  element  of  a vector  space  [this  vector  space 
is  the  tangent  space  to  the  manifold  at  the  point  p ] even  though 
q itself  is  not  an  element  of  a vector  space. 

Equation  (2.15)  is  really  a change  of  coordinates  in  the  velocity 
space  Hence,  0/  is  invertible;  calling  the  inverse  3 yields 


g = Pq 


(2.16) 


where  3 = a 1 . If  the  components  of  the  velocity  are  time 
derivatives  of  coordinates  [ e.  g.  , when  a is  the  identity  ],  then 
g is  called  a holonomic  velocity  [ strictly,  the  components  of  g are 
holonomic  velocity  components  ] ; holonomic  components  are  also 
called  integrable  components.  If  the  components  of  the  velocity 
are  not  time  derivatives  of  coordinates  [ i.  e.  , the  components  are 
not  integrable  ] , then  they  are  called  nonholonomic  velocity 
components.  Interesting  discussions  of  nonholonomic  dynamics 
equations  (also  known  as  Boltzmann-Hamel  equations  and  as 
Lagrange  equations  in  quasi-coordinates)  can  be  found  in  [l4]-[l5]. 

After  a basis  is  introduced  p , v , and  y become  n x n matrices. 
M-  and  v depend  nonlinearly  on  q:  p = p (q)  and  v=  v(q);  y depends 
nonlinearly  on  q and  linearly  on  g or  G:  y - y(q,g)  or  y = y (q,G). 

D.  Plant 


Using  the  velocity  formulation  for  the  dynamics  equation,  the 
plant  equations  can  now  be  written  in  the  form  [ see  Fig.  1 ] 
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Controller  Dynamics  Kinematics 


a)  Momentum  Formulation 


Controller  Dynamics  Kinematics 


e = 0 (e,  q,  t) 

K 

g = Vg  + VK 

g 

q = ag 

K = e,  q,  t) 

g = g 

q = q 

I 


b)  Velocity  Formulation 


Figure  1.  Block  Diagram  of  Spacecraft  Attitude 
Controller,  Dynamics,  and  Kinematics 
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These  equations  have  the  general  form  (see  Fig.  2) 

x = ^ x + u 

y = Cx 


(2. 17a) 


(2. 17b) 


(2.18a) 

(2.18b) 


Alternatively,  using  the  momentum  formulation  for  the  dynamics 
equation,  the  plant  equations  can  be  written  as 


(2.19a) 


(2.19b) 

(2.19c) 


q 


Controller 


Plant 


e = 6 (e,  q,  t) 

• 

K 

x = lAk  + 15K 

K = t (e,  q,  t) 

q = C x 

f 


a)  Momentum  Formulation 


Controller  Plant 


e = 0(e,q,t) 

K t 

x = + ^BK 

K = f (e,q,t) 

q = Cx 

l 


b)  Velocity  Formulation 


Figure  2.  Block  Diagram  of  Spacecraft  Attitude 
Controller  and  Plant 
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o'* 1 


which  have  the  form 


x = Jx  + ^Bu 


(2.  20a) 


y = 


Cx 


(2.  20b) 


A straightforward  calculation  shows  that  equations  (2-20)  can  be 
obtained  from  (2.18)  as  follows 


where 


x = T*  x 

J = (qM  + 
q?  = 

C = C ^"P"1 


(2.  21a) 
(2.  21b) 

(2.  21c) 

(2.  21d) 

(2.  21e) 


Equations  (2.  21)  have  the  same  form  as  the  transformation  equations 
for  linear  time-variable  systems  [5],  but  here  the  transformation 
matrix  T*  depends  on  the  state  variable  [ actually,  ^P  only  depends 
on  the  output  variable,  q ] . 

Since  o and  v are  n x n nonsingular  matrices,  all  the 
coefficient  matrices  ( , C » and  , Q , C ) have  rank  n. 

In  fact,  these  coefficient  matrices  are  analytic  functions  of  the  state 
variables. 

The  next  section  discusses  equations  (2.17)  to  (2.20)  from  a 
system  theory  point  of  view. 
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III.  APPLICABLE  SYSTEM  THEORY 


The  previous  section  provides  the  working  knowledge  necessary 
to  relate  system  theory  concepts  to  spacecraft  attitude  dynamics 
problems.  In  particular,  the  concepts  of  controllability,  observability, 
and  other  algebraic  system  theory  concepts  can  now  be  meaningfully 
applied.  In  addition,  the  equations  of  motion  are  now  in  a form  such 
that  the  terms  linear  in  the  state  variables  can  be  conveniently 
separated  from  the  terms  which  are  quadratic  and  higher. 

Equations  (2.18a)  and  (2.  20a)  show  that  the  general  spacecraft 
attitude  dynamics  plant  has  a system  state  equation  of  the  following 
specialized  form 

x = J(x,t)x  + cg(x,t)u  (3.1) 

where  and  are  bounded  for  bounded  x.  Furthermore,  ^4  and  T? 
are  actually  analytic  in  x [except  at  various  discrete  times  when  the 
order  or  the  structure  of  the  equations  might  change  ].  Equation  (3.1) 
is  linear  in  the  control  u (=K);  this  can  be  emphasized  by  writing 

x = F(x,t)  + q?(x,t)u  (3.2) 


where 


F(x,t)  = ^(x,t)x 


(3.3) 


with 


F(0,  t)  = 0 


(3.4) 


Equations  of  the  form  of  (3.2)  have  been  investigated  in  [ 16]  - [ 20] 
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Equations  of  the  form  in  (3. 1 ) have  been  studied  by  Davison  and 
Kunze  [21 J in  terms  of  the  related  time-varying  linear  system 
defined  by 

x = ^(z(t),  t)  x + Q(z(t),  t)  u (3.5) 

where  z(t)  is  a prescribed  function  from  the  set  of  possible  solutions 
to  Equations  (3.  3)  and  (3.  5). 

A.  Controllability 

Theorem  1 of  [ 21  ] deals  with  global  controllability  and  requires 
that  the  elements  of  A and  Q be  bounded  for  all  possible  z(t).  But  the 
^ of  equation  (2.  20a)  is  obviously  bounded  since  it  is  a constant.  The 
of  equation  (2. 18a)  is  also  bounded  since  v is  bounded  in  multi-body 
dynamics  problems;  in  fact,  the  elements  of  v are  bounded  functions 
of  coordinates  which  are  themselves  bounded.  However,  neither  the 
A of  equation  (2. 18a)  nor  the  of  (2.  20a)  are  bounded  since  each 
depends  on  y which  is  linear  in  g or  G (and  nonlinear  in  q). 

Theorem  2 of  [21  ] deals  with  local  controllability  which  is 
obtained  when  A and  are  bounded  over  some  bounded  set  of  solutions 
z(t).  Controllability  of  the  nonlinear  system  (3.1  ) is  reduced  to 
controllability  of  the  related  linear  time-varying  system  (3.  5).  But 
controllability  of  linear  time-varying  systems  is  easily  checked  via  the 
Silverman-Meadows  controllability  matrices  [ 22  ] - [24] 

<£.  = [ (3.6a) 


where 
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33  ^ = 1? 

33^  = 33  - -JT? 

[m+  1 ] _ /d  , \ ^.[m]  _ 

33  ~ I - -3  ) 33  t ^ 0,1,  . 


(3.6b) 


(3.  6c) 


. .2n-2  (3.  6d) 


For  a system  of  order  2n: 


® S ?2n 


If  i is  the  lowest  index  for  which  rank  33  ^ = rank  33 . ^ = rank  33  , 
then  i is  called  the  controllability  index  [24]  . If  rank  33  ^ = 2n 
(for  a system  of  order  2n),  then  the  system  is  controllable. 

Velocity  Formulation 

In  the  velocity  formulation,  Equations  (2.17a)  and  (2.18a)  yield 


(3.7) 


Therefore,  from  equations  (2.10b)  and  (3.7) 


33[l]  = 33  = 


It  follows  that 


-y  *v  - vy  I - y*v  _ vy 


ffv  I - a 


(3.8) 


33  33 


I O -a 

L n 


(3.9) 


Since  v and  a are  nonsingular,  it  follows  that  Rz  has  rank  2n, 
and  therefore  the  related  linear  system  is  controllable,  and  the 
controllability  index  is  2. 


Momentum  Formulation 


In  the  momentum  formulation,  equations(2. 19a)  and  (2.  20a)  yield 


A = 

‘ y 

o 

n 

9 

V = 

m 1 

n 

a 

O 

n 

O 

n 

Since  R is  constant,  it  follows  that 


Therefore 


[ R R [1  ]]  = 


(3.10) 


(3.11) 


(3. 12) 


Since  a is  nonsingular,  it  follows  that  Rz  has  rank  2n,  and  the 
momentum  formulation  is  also  controllable  (as  expected).  Note  that 
R 2 and  Rz  are  related  by 


*2  * *<52 

as  they  should  [23  ] ; here  R is  given  by  equation  (2.  21e). 


(3.13) 
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B. 


Observability 


Observability  of  the  related  time-varying  linear  system  can  be 
checked  similarly  via  the  Silverman-Meadows  observability  matrices 
[22]  - [24] 


where 


£j 


c[0] 

r[1] 


(3. 14a) 


c[0]  = o 

etl]  = c+cj  ==»  e 1 *=  c*  + jt-  0* 


(3. 14b) 
(3. 14c) 


m = 0,  1,  . . . , 2n-2 


(3. 14d) 


For  a system  of  order  2n: 


0 


(3. 14e) 


If  i is  the  lowest  index  for  which  rank  C.  = rank  0.  = rank  C , 

-l  —l+l  — 

then  i is  called  the  observability  index  [24],  If  rank  = 2n 
(for  a system  of  order  2n),  then  the  system  is  observable. 

For  both  the  velocity  and  the  momentum  formulation,  the 
output  is  q and  therefore 

('  = [ On  ln  ] = C (3. 15) 
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This  matrix  is  constant  and  of  rank  n.  C^A  and  O^A  are  also  of 
rank  n,  and  the  observability  index  is  2.  For  the  velocity  formulation 


• 

c 

c 

m 

n u J 

CuA 

. 

-2 


For  the  momentum  formulation 


£2  ■ 


n 


n 


n 


(3.16) 


■ 

• « 

- 

- 



O 

1 

c 

e 

n 

n 

cCl] 

c-J 

5 

O 

. 

n 

(3.17) 


Evidently,  both  and  have  rank  2n,  and  they  are  related  by 


—2  = £2  * 


-1 


(3.18) 


as  they  should  [23]. 

C.  Covariance,  Contravariance,  and  Duality 

The  duality  between  controllability  and  observability  is  similar 
to  the  duality  between  covariant  and  contravariant  tensors.  This 
conclusion  has  been  stated  recently  by  Hermann  and  Krener  in  the 
wider  context  of  nonlinear  differentiable  systems  [25]  . This 
duality  will  be  illustrated  here  only  for  time-variable  linear  systems. 
Consider  the  linear  time-varying  system  given  by 


x = J(t)  x + R(t)  u 
y = C(t)  x 


(3.19a) 

(3.19b) 
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Let  § (t,tQ)  be  the  transition  matrix  so  that  given  the  state  at  tQ, 
x(tg),  the  state  at  t is  given  as 


x (t ) = $ (t,  tQ) 


x(t0)  + J $ (tQ*  T)  15(T)u(T)dT 


(3.  20) 


The  system  is  controllable  if  and  only  if  the  rows  of  $ (tQ,  . ) T?(.  ) 
are  linearly  independent  functions.  This  condition  is  equivalent 
to  the  nonsingularity  of  the  controllability  Gramian  matrix  [24] 

h 

M(t0*V  = / Mt0»T)  <:R(T)  <£ (t)  * t(t0,T),lT  (3.21) 

t0 

The  relationship  between  covariance  and  controllability  is  now 
brought  out  by  the  fact  that  M(tQ  , t^)  satisfies  the  linear  matrix 
differential  equation  [5] 

-a^Mft.tj)  = ^(tJMft.tj)  + M(t,tj)  _/(t)  - q?(t)  ^(t)  (3.22) 


Ignoring  'I?,  this  equation  is  similar  to  equation  (2.  10a)  with  M and 
A taking  the  roles  of  g.  and  y , respectively. 

The  system  is  observable  if  and  only  if  the  columns  of 
C(*  ) M*  * ) are  linearly  independent  functions.  This  condition 

is  equivalent  to  the  nonsingularity  of  the  observability  Gramian 
matrix  [24] 


N(t 


O'V  = J $t(T*  t0)  CV)  C(T)*(T.t0) 


dT 


(3.  23) 
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The  relationship  between  contravariance  and  observability  is  now 
brought  out  by  the  fact  that  Nftg.tj)  satisfies  the  linear  matrix 
equation  [5] 

3fN  (t.tj)  = -Mtmt.tj)  - Nft.tjJjKt)  - cV)  C(t)  (3.24) 

Ignoring  C > this  equation  is  similar  to  equation  (2.10b)  with  v and 
7 replacing  N and  <A,  respectively. 

D.  Transformation  of  Input 

Despite  the  duality  between  controllability  and  observability,  it 
appears  that  there  is  a lack  symmetry  because  y enters  in  ^ and 
(equations  (3.9)  and  (3.12),  respectively),  whereas  y does  not 
enter  into  and  0^  (equations  (3.16)  and  (3.17),  respectively. 

This  dependence  on  y can  be  formally  eliminated  by  redefining  the 
input  as  follows : 

k'  = k - ytg  = v (K  - p ylg  ) (3.  25a) 


K'  = K + yG 


(3.  25b) 


Equation  (2.17a)  then  becomes 


Jn  k' 


(3.26) 


In  effect,  the  y in  has  been  set  to  zero,  and  hence  the  y in  eg.,  is 
also  zero.  Of  course,  the  rank  of  ^ is  still  2n,  independent  of  y . 
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Similarly,  equation  (2.19a)  becomes 
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(3.27) 

_ q _ 
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O 
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O 

_ n J 

Again,  in  effect  the  y in  J has  been  set  to  zero;  this  then  results  in 
the  y in  being  set  to  zero,  with  rank  of  remaining  2n. 

Redefining  the  input  variable  has  another  important  effect; 

J and  ^4  are  now  bounded  (since  the  unboundedness  was  due  to 
7 ).  Consequently,  by  Theorem  1 of  [ 21  ] the  nonlinear  systems  (3.26) 
and  (3.  27)  are  globally  controllable. 


E.  Perturbed  Linearized  System 


Equations  (2.17a)  and  (2.19a)  are  actually  even  controllable 
when  they  are  linearized.  Linearization  involves  dropping  the 
term  in  y and  evaluating  O'  , a , and  v at  q = 0 . More 
generally,  the  equations  can  be  written  as  these  linearized  terms, 
plus  higher  order  "perturbations".  Thus,  equation  (2-17a)  becomes 
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(3.  28) 


This  equation  has  the  general  form 


x = A X + -Bu  + h (x,  u) 


(3.29) 


with  ^4  and  hB  being  constants;  h can  be  further  expanded  as 
h (x,  u)  = h^(x)  + h^xju 

where  h2  is  quadratic  (or  higher)  in  x and  hj  is  linear  (or  higher)  in  x. 
Similarly,  equation  (2.19a)  becomes 


G O O G IK 

n n n 


q a(0)  On  q 


O G 
n 


0/  - 0/  (0)  On  q 


This  equation  has  the  general  form 


(3.  30) 


^x  + q?u  + h (x) 


(3.31) 


with  ^4  and  1?  being  constants,  and  with  h being  quadratic  in  x. 

It  is  evident  that  the  linearized  parts  of  (3.  29)  and  (3.  31)  are 
controllable.  Systems  of  this  type  are  studied  in  [26]  - [27].  A 
slightly  more  general  perturbed  system  is  studied  in  [28].  Specific 
studies  of  quadratic  systems  are  in  [29]  - [33]. 


F. 


Nonlinear  Algebraic  System  Theory 


Loosely  put,  algebraic  system  theory  is  the  theory  of  systems 
in  terms  of  the  matrices  ^B,  C rather  than  in  terms  of  the 
solutions  of  the  system  equations.  In  particular,  the  state  transition 
matrix  does  not  appear  explicitly  in  this  theory.  Much  of  the  initial 
work  on  algebraic  system  theory  is  due  to  Kalman  and  coworkers 
[34]  - [35]  . The  work  of  Silverman  and  Meadows  [22]  is  also  part 
of  algebraic  system  theory  (for  time-variable  linear  systems)  since 
it  deals  with  controllability  and  observability  in  terms  of  _A  , T?  , C 
rather  than  in  terms  of  the  state  transition  matrix. 

Algebraic  system  theory  is  presently  still  only  a linear  theory, 
but  it  is  conjectured  that  many  of  the  results  which  depend  only  on 
^ , C can  be  extended  to  nonlinear  systems  like  equation 
(3.  3).  For  example,  the  theory  of  "linear  algebra  with  continuously 
parametrized  elements"  in  [36]  - [37]  can  probably  be  extended  to 
a theory  of  "linear  algebra  with  continuous  state-dependent  elements". 
It  is  conjectured  that  such  a theory  could  prove  that  the  of 
equation  (3.  7)  can  be  transformed  so  that  y becomes  zero: 
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(3.  32) 


where  or  is  a nonsingular  n x n matrix,  so  that  the  rank  of  the  new 
is  still  n.  The  ^ of  equation  (3. 10)  can  be  transformed 
similarly: 


O 


n 


O 

n 


(3.  33) 
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where  a is  also  a nonsingular  n x n matrix  so  that  the  rank  of 
of  the  new  A is  still  n. 

It  is  not  really  clear  whether  or  not  a change  of  variables 
can  be  found  so  that  y is  transformed  to  zero.  Even  if  such  a 
transformation  does  exists,  it  may  not  be  directly  useable  since  it 
may  depend  on  knowledge  of  the  solution  of  the  differential  equations. 
However,  it  follows  from  general  tensorial  considerations  that  the 
form  of  y can  be  simplified. 

Recall  that  y depends  nonlinearly  on  q and  linearly  on  g: 
y ~ y (q>  g)-  If  the  matrix  elements  of  y are  linearly  expanded 
in  terms  of  the  elements  of  g,  then  the  expansion  coefficients 
are  the  components  of  the  affine  connection  (for  the  chosen  basis)[l3]; 
if  the  elements  of  g are  holonomic  velocity  components  then  these 
coefficients  reduce  to  the  Christoffel  symbols  of  the  second  kind 
for  the  metric  matrix  |X  [ll]  . But  it  is  known  that  transforming 
to  geodesic  coordinates  (or  normal  coordinates)  [ll]  at  the  point  p 
yields  Christoffel  symbols  which  vanish  at  the  point  p . If  q 
represents  the  coordinates  of  the  point  p , then  it  follows  that 
y ( 0,  g ) = 0;  thus,  y is  linear  (or  higher)  in  q,  and  as  a result, 
y is  quadratic  (or  higher)  in  the  state. 

Another  interesting  fact  about  y is  that  it  can  be  transformed 
so  as  to  become  skew- symmetric.  This  is  evident  on  examining 
equations  (2.  la)  and  (2.  3a),  and  then  transforming  these  equations 
to  G = yG  + K and  g = - y*g  + k , respectively.  If 
these  transformations  are  made  such  that  G = g (i.  e.  , momentum 
and  velocity  are  "normalized"  by  the  square-root  of  the  mass)  then 
it  follows  that  y — - y , which  is  the  condition  for  skew- symmetry. 

In  order  to  be  sufficiently  general,  algebraic  system  theory 
must  be  formulated  in  terms  of  manifolds  rather  than  vector  spaces 
because  the  configuration  space  for  a dynamical  system  is  a manifold, 

Ak  , and  the  state  space  is  the  tangent  bundle,  <0T(<J4  )*  of  this  mani- 
fold. An  interesting  introduction  to  differentiable  manifolds  in  nonlinear 
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system  theory  is  the  review  by  Brockett  [38]  . More  advanced 
material  is  in  [25]  and  [39]  - [40]  . 

The  next  section  casts  multi-body  dynamics  in  a nonlinear 
algebraic  system  theory  setting.  The  multi-body  system  is  considered 
to  consist  of  N interconnected  subsystems,  and  the  nonlinear 
equations  of  motion  for  the  system  are  algebraically  related  to  the 
nonlinear  equations  of  motion  for  the  N subsystems.  The  inter- 
connection constraints  among  the  N subsystems  are  linear. 

Therefore,  these  constraints  decompose  the  (linear)  velocity  and 
momentum  spaces  into  "free"  and  "constrained"  parts;  these  are, 
in  effect,  the  controllable  and  uncontrollable  subspaces.  "Short 
exact  sequences"  and  "algebraic  network  diagrams"  are  used  to 
relate  the  various  subspaces. 
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IV.  MULTI-BODY  DYNAMICS 


| 

i 


Consider  a (dynamical)  system  which  can  be  conveniently 
considered  as  a collection  of  N interconnected  subsystems.  For 
example,  in  classical  dynamics  it  is  well  known  [15  ] that  any 
mechanical  system  can  be  considered  as  a collection  of  N 
particles.  Each  particle  has  three  (translational)  degrees  of 
freedom,  so  that  the  system  has  a total  of  3N  degrees  of  free- 
dom if  the  interconnections  do  not  impose  any  constraints.  If 
the  interconnections  impose  m constraints,  then  there  are 
n = 3N  - m degrees  of  freedom;  the  resulting  dynamics  plant  is 
of  order  2n,  having  n free  velocity  components  and  n free  position 
coordinates. 

Often  a mechanical  system  can  be  considered  as  a collection 
of  N rigid  bodies  [l]  - [4];  each  rigid  body  has  six  degrees  of 
freedom  (three  translational  and  three  rotational),  so  that  the 
system  has  6N  degrees  of  freedom  before  the  application  of 
constraints.  After  the  incorporation  of  m constraints,  there  are 
n = 6N  - m degrees  of  freedom;  the  resulting  plant  is,  of  course, 
again  of  order  2n.  More  generally,  a mechanical  system  can  be 
considered  as  a collection  of  N deformable  bodies,  with  Body  i 
having  n^^  s 1 degrees  of  freedom. 


A.  Primitive  Composite  System 


Now  consider  a collection  of  N bodies,  with  Body  i having  n^ 
degrees  of  freedom.  The  momentum  formulations  dynamics  equations 
for  Body  i are 


'(i) 


— y...  G...  + K... 

(l)  (l)  (l) 


(4.1a) 


’(i) 


= V(i)  G(i) 


(4.  lb) 


Here  K^,  ^(i)’  anc^  ^(i)  are  co^umn  matrices  with  n^  elements, 

and  y...  and  v...  are  n...  x n...  matrices.  The  velocity 
(i)  (i)  (i)  (i)  1 

formulation  dynamics  equations  for  Body  i are 


y(i)  g(D  + 


V(i)  Kd) 


(4.2a) 


(4.  2b) 


K^),  , and  g^  are  called  the  Body  i force,  momentum,  and 

velocity,  respectively;  v and  are  the  Body  i inverse  mass 

matrix  and  Cartan  matrix,  respectively. 

Ignoring  the  constraints,  the  system  force,  momentum  and 
velocity  are  given  in  terms  of  direct  sums  of  the  corresponding 
subsystem  quantities,  as  follows: 


, G = 


(4.  3a,  b,  c) 


Similarly,  the  system  inverse  mass  matrix  and  Cartan  matrix  are 
given  in  terms  of  direct  sums  as  follows 


(4.  3d) 
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y = 


^(N) 


(4.3e) 


Note  that  v and  y are  "block-diagonal",  with  the  off-diagonal 
blocks  being  zero. 

When  the  constraints  are  ignored,  the  momentum  formulation 
system  equations  of  motion  are  given  by 

G = y G + K (4. 4a) 


g = vG 


(4.  4b) 


and  the  velocity  formulation  system  equations  are 

g = -yt  g + VK 

g = g 


(4.5a) 
(4.  5b) 


These  equations  are,  of  course,  exactly  the  same  as  equations  (2.1) 
and  (2.  3). 

The  system  mass  matrix  is  the  direct  sum  of  the  subsystem 
mass  matrices: 


(4.6) 


-35- 


M<  is  the  inverse  of  v , and  hence 


G = Mg 


(4.7) 


as  in  equation  (2.  2). 

Let  denote  the  kinetic  energy  of  Subsystem  i.  Then 

and  satisfy 

T(i)  = 7 G(i)g(i)  = 7g(i)^(i)g(i)  = rG(i)v(i)G(i)  (4*8a) 


T(i)  " K(i)g(i) 


(4.8b) 


The  kinetic  energy  of  the  system,  with  the  constraints  ignored,  is 
then  defined  by 

N 

t * Z T(ii  <4-  9> 


i=l 


Thus  T and  T satisfy 


rp  _ 1 „t  It 

T - 2Gg_2g  ^g 


J G*  v G 


(4.10a) 


T = Kg 


(4.10b) 


as  in  equations  (2.  6)  and  (2.  9). 

When  the  constraints  are  ignored,  the  system  is  called  the 
primitive  composite  system  since  this  system  is  in  fact  a collection 
of  N independent  and  non-interacting  subsystems.  The  primitive 
composite  system  is  also  called  the  "torn  system",  after  Kron  [ 4l] . 
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The  primitive  composite  system  equations  are  converted  to  the 
actual  or  constrained  composite  system  equations  by  incorporation 
of  the  constraint  equations.  But  before  doing  this,  it  is  of  interest 
to  examine  how  the  equations  transform  under  a linear  velocity 
(or  force)  transformation.  This  velocity  transformation  can  then 
be  selected  to  impose  the  required  constraints. 

B.  Linear  Velocity  Transformation 

Consider  the  nonsingular  linear  velocity  transformation 

g = A g (4.11) 

where  A = A(q)  depends  on  the  position,  but  not  on  the  velocity. 

Let  B be  the  inverse  of  A: 

B = A'1  (4.12) 

Equation  (4.11)  can  now  be  inverted  as  follows: 

g = B g (4.13) 

The  velocity  transformation  A induces  momentum  and  force 
transformations  by  A*  as  follows: 

G = A*  G (4.14) 

and 

K = A4  K (4. 15) 

These  are  the  required  definitions  for  G and  K in  ordei  for  T and 
T to  have  the  same  form  in  terms  of  both  the  barred  and  unbarred 
quantities: 

T = j ctg  = j G*g  (4.16) 

T = Ktg  = K*g  (4.17) 
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In  order  to  get 


_ _ It  1 _t 

T - 2 g g = 2 g ^ g (4.18) 

it  is  necessary  to  define 

M-  = A1  M-  A (4.19) 

There  now  results  the  relationship 

G = M.  I (4.20) 

which  has  the  same  form  as  G = (Jt  g . 

Equations  (4.14)  to  (4.20)  also  have  their  inverse  counterparts. 

Thus , 

G = 13*5  (4.21) 

K = B*K  (4.22) 

and 

T = I G*  v G = | G vg  (4.  23) 

where 

v = B v b1  = ]I  _1  (4.24) 

Equation  (4.  20)  can  now  be  inverted  as 

g = v G (4.25) 

which  has  the  same  form  as  g = v G . 

The  barred  momentum  formulation  equation  of  motion  is  now 
given  by 
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G - y G + K (4.  26a) 

where 

y = At  y Bl  + A1  B*  (4.  26b) 

The  barred  velocity  formulation  equation  of  motion  is  now 
given  by 

i = -y*g  + k (4.27a) 

where 

k = Bk  (4.27b) 

and 

9t  = B ytA  - BA  (4.27c) 


Equation  (4.  27c)  is  equivalent  to  (4.  26b)  since 

B A = !n=^  B A = - B A (4.  28a,  b) 

Equation  (4.  27b)  also  has  the  inverse  form 

k=  Ak  (4.  29) 

which  is  analogous  to  g = A g . From  K = A^  and  K = g.  k there 
now  follows 

K = jik  k = v K (4.  30a, b) 

It  is  interesting  to  note  that  all  the  transformation  equations 
have  the  same  form  as  would  be  obtained  if  A (and  B)  were  time- 
dependent  rather  than  position  dependent.  Thus,  the  transformation 
theory  of  linear  systems  carries  over  more  or  less  directly  to 
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I 

I 


I 


i 


I ' 


I 

I 


I 


nonlinear  systems  with  linear  velocity  (or  momentum)  transformations. 
It  should  also  be  noted  that  the  coordinate  transformation  (implied 
by  the  above  velocity  transformation)  is  not  shown;  this  coordinate 
transformation  is  in  general  nonlinear  since  the  configuration 
manifold  is  a nonlinear  space.  But  the  associated  velocity  and 
momentum  transformations  are  linear  since  the  tangent  and  cotangent 
spaces  at  a point  p are  linear  vector  spaces. 


C.  Partitions  and  Decompositions 


Generally  A and  g are  chosen  so  that  the  resulting  g 
can  be  partitioned  into  a part  g^  which  is  "free"  or  "unknown"  and 
a part  g2  which  is  "constrained"  or  "known"; 


g 


(4.  31) 


This  partitioning  of  g induces  a corresponding  partitioning  of  A: 

A = [ Aj  A2  ] (4.  32) 


so  that  the  equation  g = A g becomes 


g = C Aj  A2  ] 


§1 

h 


= Aih  + A2g2  <4-  33> 


Now  suppose  that  g^  has  n^  elements,  g2  has  n2  elements,  and  g 
has  n = n.  + n_  elements.  Then  A.  is  an  n X n.  matrix  with 

1 d i i 

rank  m,  for  i = 1,  2. 
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The  partitioning  of  g also  induces  a corresponding 
partitioning  of  B: 

Tb,  1 


so  that  the  equation  g = B g becomes 

r *ii  r®i * 


LM  LBz  gJ 

Since  A and  B are  inverses  of  each  other  , it  follows  that 


(4.34) 


(4.35) 


A B = A1B1  + A2  B2  = 1 


where  1 is  the  n X n identity.  Also 


B1  A1  B1  A2  *1  ° 


B2  A1  B2  A2  ° lZ 


(4.36) 


(4.37) 


where  1.  is  the  m x m identity.  Thus  B^  is  a left  inverse  of 

A^,  and  A^  is  a right  inverse  of  B^,  for  i = 1,  2. 

Since  g^  is  free  whereas  g ^ is  constrained,  equation  (4.  33)  is 

a unique  decomposition  of  g into  its  free  and  constrained  parts. 

1 2 

Define  g and  g by 


g = A.  g.  , 1 = l>  2 


(4.  38) 


Then  equation  (4.  33)  becomes 


1 X 2 

g - g + g 


(4.  39) 


Next,  note  that 


Bi  g = Big  = gA  » i = 1,2 


(4.  40) 


The  momentum  can  be  partitioned  and  decomposed  analogously. 
From  G = AG  it  is  evident  that  G is  partitioned  into  G^  and  G ^ 
as  follows: 


G = 


f 

G1 

aJ  g 

- 

t 

G2 

A2g 

_ _ 

_ 

(4.  41) 


From  G = B*G  there  follows 


G = B|  G,  + B*  G2 


(4.42) 


Thus,  G is  decomposed  as 


s 2 

G = G + G 


(4.43) 


where 


G1  = B*  G.  , i a 1,2 

li 


(4.44) 


-42- 


Next,  note  that 


At  G = At  G1  = G.  , i = 1, 2 
1 l l ’ 


(4.45) 


The  foice  can  be  partitioned  and  decomposed  analogously. 
— t — — i- 

From  K = A K it  is  evident  that  K is  partitioned  into  and 

as  follows: 


— — 

— — 

t 

K1 

AJ  K 

1 

(VJ 

Itt 

1 

a'k 

(4.46) 


K.  is  the  "free"  or  "control"  force,  whereas  K~  is  the  constraint 

£ ^ 

force.  From  K = B K there  follows 

K = bJ  Kl  + B*  K (4.47) 

Thus,  K is  decomposed  as 

K = K1  + K2  (4.48) 

where 

K1  = b‘  K.  , i = 1 , 2 (4.  49) 

Next,  note  that 

A*  K = A*K*  = K.  , i = l,2  (4.50) 

Note  that  the  sequence  of  transformations  could  have  been  started  with 
selecting  A and  K so  that  the  resulting  and  K2  have  the  above 
properties.  The  velocity  transformation  which  specifies  g^  and 
would  then  be  an  induced  transformation  rather  than  the  primitive 
transformation. 
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The  kinetic  energy  T can  now  be  expressed  as 


T 


1 

2 


-t 

G g 


1 ir1 

2 (G 


A i f"  “ “ \ 

g + G g ) 


(4.  51a) 


= I <G'l  6!  + G282> 


(4.  51b) 


and  the  time  derivative  of  the  kinetic  energy  becomes 


T 


1 

g 


+ 


2 

g 


(4.  52a) 


= k'Ij  + K*  g2 


(4.  52b) 


Note  that  the  "constrained"  (or  "prescribed")  velocity  g^  contributes 
to  T and  T,  unless  g^  is  zero.  From  (4.  52b)  it  is  evident  that  if 
g^  is  non-zero  then  the  constraint  force  performs  work  and  thus 
affects  the  system's  energy. 

D.  Induced  Subspaces  and  Short  Exact  Sequences 
Now  introduce  the  notation 

tt1  = A.B.  , i = 1,2  (4.53) 

From  Equation  (4.  36)  it  follows  that 

tt1  + n2  = 1 (4.54) 

Since  Ek  is  a left  inverse  of  it  follows  that  tt1  is  idempotent: 

TT1  TT1  = TT1  , i = 1,  2 (4.  55) 
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Also 


TT 


1 


2 

TT 


O = TT 


2 


1 

TT 


(4.  56) 


It  follows  that  tt  and  tt  are  orthogonal  projectors  [42]  . These 
projectors  decompose  the  velocity  space  qy  into  direct  sums  Qy 
and  q/ 2 with  g € qy  and  g*  C qy*  : 

qy  = qy1  © <v2  (4. 57) 

where  © denotes  a direct  sum  of  vector  spaces.  Similarly,  let 
g C q/  and  g^  C qy.  ; then  qy  is  the  direct  sum  of  q^ , and  qy^  : 

qy=  qyx  © qy2  (4.58) 


Since  g and  g are  related  by  a nonsingular  linear  transformation,  it 
follows  that  qy  and  qi  are  isomorphic: 


qy  ~ qy 


(4.59) 


.t 


The  projectors  n1  decompose  the  dual  space  Q = 

12  • 
direct  sums  Q and  Q with  G f (j  and  G*  C Q1  : 


qy*  into 


9 = g1  © g2  (4.60) 

Similarly  9 = qy  is  decomposed  into  direct  sums  9^  and  92  with 
G € 9 and  G^  € : 


g = 9j  © 92  (4.6i) 
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Again,  Q and  Q are  isomorphic: 


9^9  (4.62) 


The  partitions  of  A and  B can  be  used  to  form  four  short 
exact  sequences  [43]  . Two  of  these  sequences  deal  with 
and  <^2  : 


A1 

B2 

q^l  - 

► a 

^ ^2 

(4.  63a) 

A2 

B1 

^2  - 

(4.63b) 

Two  more  sequences 

deal  with 

9 ^ and 

9 2 : 

BI 

A1 

2 

" 

9 

(4.  64a) 

B2 

A1 

A1 

- 

9 

% 

(4.  64b) 

Equations  (4.  63a)  and  (4.  64b)  can  be  combined  as  follows 


q/ 


^i 


r 

9, 


V. 


22 


(4.65) 


9 


9, 
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Similarly,  equations  (4.63b)  and  (4.64a)  can  be  combined  as  follows: 


B, 


22  ; 


- <\) 


- q/i 


T 

9o 


9 


A 


9, 


li 


(4.66) 


These  last  two  diagrams  look  like  the  "algebraic  network  diagrams" 
of  Roth  [44]  - [45]  and  Branin  [46]  - [47]  . These  diagrams  can 
be  filled  in  as  shown  via  dashed  lines.  The  notation  for  lL..  and 

i’ 

, for  i = 1,  2,  shown  in  these  diagrams  is  obtained  from  the 
partitioned  forms  of  n and  v . Thus 


M.  = A n A = 


Al  M A: 


A2  A]_ 


A1  “ A2 


A2  » A2 


“11 

^12 

S21 

^22 

(4.67) 


Similarly, 


v = B v B = 


B1V 

B1 

B.  v b! 

1 l 

B2v 

BI 

B2  v B 

11 

V12 

21 

V22 
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(4.68) 





Since  p and  v are  positive  definite  symmetric,  and  and 
have  full  rank,  it  follows  that  and  <re  also  positive 

definite  symmetric,  and,  therefore,  invertible. 

From  (4.65),  for  example,  it  should  not  be  concluded  that 
is  equal  to  ; this  is  so  only  if  g2  is  zero.  In  general, 

from  G = pg.or  from  (4.  65)  plus  (4.66)  there  follows 


The  inverse  relationship, 


G1  = 

“u  §i 

+ “l2  ®2 

(4.  69a) 

G2  = 

“21  8l 

+ “22^2 

(4.  69b) 

g = 

v G , 

yields 

6l  = 

*11  G1 

+ V 12  G2 

(4.  70a) 

h = 

V21  G1 

+ ~22  G2 

(4.  70b) 

E.  Constrained  Motion  Equations 

Using  equation  (4.47)  for  K,  the  momentum  formulation  dynamics 
equation  can  be  written  as 

G = y G + b|  Kj  + B*  K2  (4.  71) 

Here  K.,  is  the  unknown  constraint  force  which  can  be  determined  via 
use  of  the  constraint  equation  B2g  = g2  . Writing  this  constraint 
equation  in  terms  of  momentum  yields 

B2  v G = g2  (4.  72) 
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Differentiating  this  equation  with  respect  to  time  and  then  combining 
with  equation  (4.  71)  yields 


1 

n 

-B2 

G 

y G + bJ  k1 

B2V 

O 

*2 

h - B2  V G - B2  * G 





__  _ 

_ 

(4.73) 


Solving  this  for  G yields  G with  eliminated  (solving  the  equation, 
of  course,  also  yields  K^).  However,  the  resulting  equation  for  G 
is  not  a minimal  dimension  realization  of  the  dynamics,  and  hence 
the  system  is  not  controllable. 

To  get  a minimal  dimension  realization,  left-multiply  equation 
(4.  71)  by  a|  . Using 


Ai B! 


= l 


l ’ 


A*  B*' 
A1  B2 


= O 


(4.  74a, b) 


then  yields 


Aj  G = AjyG  + 


(4.75) 


Thus,  K 2 has  been  eliminated.  Differentiating  Gj  = A^  G with 
respect  to  time  yields 


Gj  = a\g  + AjG 


(4.76) 


Combining  equations  (4.75)  and  (4.76)  yields 
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4 


gi  = (A; 

y + 

a\  ) G 

+ 

K1 

(4.  77a) 

, . t 

• t . t 

= (A2 

y + 

Aj  ) B 

G 

+ K1 

(4.  77b) 

tion  G involves  both 

G^  and 

G2 

• Cl  J i S y 

of  course, 

the  momentum  state  variable  for  the  reduced  (and  controllable) 
system.  G£  can  be  expressed  in  terms  of  and  the  prescribed 
velocity  as  follows.  From  (4.69a)  is  given  as 


gl  ^11  * G1  " ^12  g2  * 


(4.78) 


C2  is  now  obtained  by  substituting  this  expression  into  (4.  69b).  The 
result  is 


G2  " ^21  ^11  G1  + ( ^22  ' ^21  ^11  ^12  * g2 


(4.79) 


If  preferred,  the  use  of  can  be  replaced  with  the  use  of 

as  follows 


IL-^  = V — V \>  ~i  \1 

HI  11  12  22  21 


(4.80) 


--1 


Note  that  the  coefficient  of  g ^ in  equation  (4.  79)  is  just  ^2  : 


V22  = ^22  ' ^1  *11  ^12 


(4.  81) 


Thus,  G2  can  also  be  expressed  as 

°2  m >*21  < ”U  - "l2  722  V21  > S1  + “2I  *2  <4'82» 
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The  constrained  motion  equations  have  been  given  only  for  the 
momentum  formulation,  but  the  velocity  formulation  equations  are 
very  similar  [3].  In  general , there  is  no  significant  conceptual 
difference  between  the  momentum  formulation  and  the  velocity 
formulation,  and  between  carrying  the  constraints  as  a side 
condition  (as  in  equation  (4.  73))  and  eliminating  the  constraints 
(as  in  equation  (4.  77)).  Whatever  approach  is  selected,  there  is 
the  choice  of  (effectively)  inverting  or  v , both  of  which 

are  positive  definite  and  symmetric.  However,  in  particular 
multi-body  systems,  one  of  these  four  approaches  may  have 
computational  advantages  over  the  others  [3]. 

The  search  for  various  alternative  formulations  of  multi- 
body dynamics  is  analogous  to  the  search  for  various  minimal 
and  "canonical"  realizations  in  system  theory.  Although  various 
realizations  may  be  equivalent,  different  points  of  view  and  different 
insights  can  be  obtained  from  different  approaches. 


V.  DISCUSSION  AND  CONCLUSION 


The  previous  three  sections  discuss  spacecraft  attitude  dynamics 
and  multi-body  dynamics  from  the  point  of  view  of  system  theory.  It 
is  hoped  that  this  discussion  is  of  interest  to  system  theorist,  who 
may  be  looking  for  a sample  problem  of  modest  complexity,  and  to 
dynamicists,  who  might  be  interested  in  the  techniques  of  system 
theory. 

It  is  expected  that  in  time  the  communication  gap  between  system 
theorists  and  dynamicists  will  narrow.  This  narrowing  is  expected 
to  occur  most  rapidly  in  the  field  of  spacecraft  attitude  dynamics 
where  the  trends  are  to  larger  size  plants  (both  mathematically  and 
physically)  and  more  sophisticated  controllers  (including  micro- 
processors in  both  the  sensors  and  control  laws).  For  a geometric 
perspective  on  the  relationship  of  control  theory  and  mechanics 
see  [48]. 

The  reader  interested  in  more  details  on  multi-body  dynamics 
is  referred  to  [l]  - [3]  and  [49]  - [56],  In  these  references  the 
momentum  formulation  equation  of  motion  has  the  general  form 

6 + X = K (5.1) 

and  the  velocity  formulation  equation  of  motion  is 

Hg  + Y = K (5.2) 

The  relationship  between  X and  Y is 


Y = X + pg 


(5.3) 


Comparing  Equations  (5. 1)  and  (5.  2)  with  Equations  (2.  la)  and  (2.  3a), 
respectively,  shows  that 

X = -yc  (5.4a) 

Y = \iytg  (5.4b) 

Equations  (5.  4)  are  useful  in  relating  the  multi-body  dynamics 
formulations  in  the  references  with  the  formulation  presented  in 
this  paper. 

Left-multiplying  Equations  (5.1)  and  (5.2)  by  v = yields 

vG  + x = k (5.5) 

and 

g + y = k (5.6) 

where 

x = vx  = - vyG  (5.  7a) 

y = v y = ytg  (5.  7b) 

k = vk  (5.7c) 

Of  course  the  variables  x and  y introduced  here  differ  from  the  x 
and  y of  Sections  II  and  III. 

Equations  (5.  2)  and  (5.  5)  are  generally  not  as  convenient  as 
Equations  (5.  6)  and  (5. 1)  , respectively,  but  they  are  given  here 
( for  completeness.  A very  practical  question  is  which  of  these 

four  formulations  is  computationally  most  convenient.  As  pointed 

■ 

out  in  Section  II,  generally  it  is  not  convenient  to  find  y since  y 
itself  is  not  required  in  the  equations  of  motion.  Similarly,  p and 
^ are  not  required  since  these  terms  always  occur  in  the  combinations 
>•  of  M-g  and  vG  . It  is  easy  to  see  that 


«• 


Y - X 


M-g  = 

fV  + 

“y8)e 

(5.8a) 

VG 

= (v'v  + 

vyc>G 

(5.8b) 

where  y has  been  decomposed  into 


and 


y = yG  + y'G  where  y^G  = 0 


y ~ y + y ' where 
g g 


y tg  = o 

g s 


(5.9a) 


(5.9b) 


It  now  follows  that  only  the  y^.  part  of  y is  required  in  the  momentum 

formulation  whereas  only  the  y part  of  y is  required  in  the 

g 

velocity  formulation. 

Examples  which  illustrate  many  of  the  ideas  and  concepts  in 
this  paper  are  contained  in  the  Appendix  of  a forthcoming  Aerospace 
Corporation  Technical  Report  (TR)  with  the  same  title  as  this  paper. 

The  work  reported  in  this  paper  is  a natural  extension  of  the 
author's  earlier  work  on  the  transformation  operator  approach  to 
multi-body  dynamics.  This  earlier  work  in  turn,  is  an  extension 
or  adaptation  of  some  basic  ideas  in  Kron's  method  of  sub~paces; 
see  [57]  and  [58,  Part  A ] . 
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APPENDIX 


The  purpose  of  this  appendix  is  to  illustrate  some  of  the  key 
ideas  and  concepts  presented  in  the  paper.  Two  examples  are 
considered.  The  first  is  a single  rigid  body  viewed  as  a system  of 

N particles.  As  a byproduct,  the  Cartan  operator  for  a rigid  body 
is  derived. 

The  second  example  is  the  interconnection  of  two  rigid  bodies. 
This  illustrates  the  fact  that  even  though  the  equations  of  motion  are 
nonlinear  differential  equations,  nevertheless,  the  interconnection 
constraint  equations  are  linear  in  velocity,  momentum,  and  force. 

In  the  following  a classical  3-dimensional  vector  (or  first 
order  tensor)  is  denoted  by  an  overbar  with  a single  arrowhead; 
a classical  3-dimensional  dyadic  (or  second  order  tensor)  is 
denoted  by  an  overbar  with  a double  arrowhead. 


A.  FROM  PARTICLES  TO  A RIGID  BODY 

A single  rigid  body  has  six  degrees  of  freedom.  The  three 
rotational  degrees  of  freedom  are  determined  by  the  rotational 
equation  of  motion 


H 

c 


(A.  la) 


where  Lc  is  the  torque  (or  moment  of  force)  about  the  center  of 

mass  c,  Hc  is  the  angular  momentum  (or  moment  of  linear 

momentum)  about  c,  and  the  dot  over  H denotes  the  inertial 

c 

time  derivative.  The  three  translational  degrees  of  freedom  are 
determined  by  the  translational  equation  of  motion 


P = F (A.  lb) 

where  F is  the  force  and  P is  the  linear  momentum.  The  momenta 
are  linearly  related  to  the  velocities  as  follows: 

Hc  = T • u>  (A.  2a) 

P = M • v 


where  I is  the  inertia  dyadic,  M is  the  mass  dyadic,  u>  is  the 
angular  velocity,  and  vc  is  the  linear  velocity  of  the  center  of 
mass  c.  I and  M are  positive  definite  symmetric  dyadics  which 
have  the  positive  definite  symmetric  inverses  J and  W , respectively: 


J = I (A.  3a) 

W = M _1  (A.  3b) 
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I 

I 


Consequently  equations  (A.  2)  have  the  inverse  relationships 

= "j  • Hc  (A.  4a) 

vc  = W • P (A.  4b) 

Equations  (A.  1)  are  the  explicit  momentum  formulation  equations  of 
the  form  G + X = K . Equations  (A.  2)  and  (A.  4)  are  the  linear 
relationships  G = M*  • g and  g = v . G,  respectively. 

The  coupled  velocity  formulation  equations  of  the  form 
g • g + Y = K are  given  by 


(A.  5a) 
(A.  5b) 

where  w is  the  skew- symmetric  dyadic  of  the  vector  uo  (so  that 

u>  • H = uj  X H ). 
c c 

The  mass  and  inverse  mass  dyadics  are  not  only  positive 
definite  symmetric,  but  they  are  also  spherical.  Thus, 


•*— ► 

M = 

M E 

(A.  6a) 

II 

W E 

(A.  6b) 

W = 

* 

M"1 

(A.  6c) 

j where  M and  W are  the  scalar  mass 

and  inverse  mass. 

respectively, 

and  E is  the  identity  dyadic. 


A.  1.  System  of  N Particles 


By  considering  the  rigid  body  to  be  a system  of  N particles, 
the  above  equations  of  motion  for  the  rigid  body  can  be  derived  from 
the  equations  of  motion  for  a particle  plus  the  interconnection 

tVi 

constraint  equations.  Let  rm  be  the  mass  of  the  i particle 

and  let  r^  be  its  position  vector  from  an  inertial  reference  origin. 

The  inertial  time  derivative  of  r.  is  the  velocity  ; thus,  v.  = r . . 

The  center  of  mass  c of  the  system  of  N particles  has  position 

vector  r where 
c 

N 

M ?c  = £ r.  (A.  7a) 

i=l  1 

where 

N 

M = £ m.  (A.  7b) 

i=l  1 

-♦ 

Let  R^c  be  the  position  vector  to  the  i particle  from  the  center 
of  mass.  Then  r^  can  be  decomposed  as  follows 


r. 

l 


= R. 

1C 


(A.  8) 


Taking  the  inertial  time  derivative  of  this  equation  yields 


v. 

l 


(A.  9) 


where  vc  = rc  is  the  velocity  of  the  center  of  mass  which 
satisfies  the  relationship 


N 

M v = Z m v.  (A.  10) 

i=l  1 1 
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The  relative  position  vector  R.  is  fixed  in  the  rigid  body. 
Hence  the  time  variation  of  R.  is  due  entirely  to  the  rotation 

1C  ’ 

with  angular  velocity  ; consequently 


R.  = ou  x R. 

1C  1C 


-R. 

1C 


U)  = 


R* 

1C 


(A.  11) 


where  R*  is  the  dyadic  transpose  of  the  dyadic  R.  Since  R is 
skew-symmetric  it  follows  that  R*  = -R. 

Equation  (A.  9)  can  now  be  written  as 


v.  = R^  • ou  + v 
1 1C  c 


(A.  12) 


From  this  there  follows  the  matrix-vector-dyadic  expression 


~t  -* 

71 

R,  E 

lc 

OU 

~t  rr 

v2 

= 

R,  E 

2c 

V 

c 

• 

~ t tr 

7n 

R'  E 

Nc 

This  equation  has  the  general  form 

g = A • g 


\ 

(A.  13) 


(A.  14) 
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where  g and  g are  column  matrices  of  vectors  and  A is  a 
matrix  of  dyadic s.  The  equation  G = A*  • G takes  the  form 


r 

H 

c 

o 

1 

R, 

2c 

• • • 

Rm 

Nc 

P 

E 

♦♦ 

E 

• • • 

E 

Pi 

P2 


(A.  15) 


i.V 

where  = rm  is  the  linear  momentum  of  the  1 particle. 

Thus,  the  equation  G = A*  • G yields  the  fundamental  results 


N 


5 = V]  R. 

c / J i 

i=l 


p. 

1C  *i 


(A.  16a) 


N 

5 ■ E \ 

i=l 


(A.  16b) 


Note  that  Hc  and  P are  the  new  momentum  variables  associated 

with  the  new  velocity  variables  u>  and  v . Similarly,  the 

. T*  t ^ 

equation  K = A • K yields  the  fundamental  results 


N 


l = r. 

C i—J  1C 


f. 

l 


(A.  17a) 


i=l 


N 


* • E h 

i=l 
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(A.  17b) 


L 


where  f.  is  the  force  on  the  i particle,  so  that  Newton's  Law 
of  Motion  for  the  particle  becomes 


m. 

1 


v. 

X 


f. 

1 


(A.  18) 


Again,  note  that  L£  and  F are  the  new  force  variables  associated 
with  the  new  velocity  variables  and  v^  . 

The  equation  G = p • g takes  the  form 


— — 

— ~ 

Pi 

E 

p2 

= 

m2  E 

• 

V2 

Pn 

mN^ 





— 

(A.  19) 


where  the  off-diagonal  elements  of  p are  zero  dyadic s O . 
The  equation  pi  = A*  • p • A now  yields  the  result 


N 

N 

E 

m.  R.  • R? 

1 1C  1C 

m.  R. 

/ ■>  1 1C 

i=l 

i=l 

N 

N 

E 

m.  R* 

1 1C 

m.E 

i=l 

i=l 

(A.  20) 


I O 

O M 

Thus,  the  equation  p = A*  • p • A yields  the  fundamental  result 


N 


I = L m.  R. 
i=l  1 1 


1C 


R* 

1C 


(A.  21) 
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Note  that  I and  M form  the  block-diagonal  dyadic  mass  operator 
associated  with  the  new  velocity  and  momentum  variables. 


A.  2 Evaluation  of  Cartan  Operator 

There  remains  the  evaluation  of  the  new  Cartan  operator  y 
from  the  expression  y = At  • y • B*  + A*  • B*  . Since 
y is  zero,  y reduces  to 

Y = A1  • B1  (A.  22) 


A is  obtained  from  Equation  (A.  13).  This  operator  is  one-to-one 
and  therefore  has  a left  inverse  (in  fact,  an  infinite  number  of  left 
inverses).  A specific  left  inverse  B is  chosen  by  letting  B be 
the  pseudo-inverse  of  A as  follows 


where 


Then 


B 


= v 


M* 


(A.  23) 


(A.  24) 


y = a1 


(A.  25) 


For  convenience,  introduce  6 


as  follows 


6 = 


(A.  26) 
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1 

4 

Performing  the  indicated  operations  yields  ! 


The  evaluation  of  D requires  several  uses  of  the  identity 


Consequently,  from  equation  (A.  11)  there  now  follows 


R.  = - R.  • u>  + u>  . R. 

1C  1C  1C 


(A.  30) 


Hence 


N 


D = Y'  m.  R.  • B.} 

A— I 1 1C  1C 

i=l 


(A.  31) 


N 


£ m.  R. 

' 1 1C 

i=l 


uu  • r : 


' + V m. 

1C  1 


(JO  • R.  • R. 

1C  1C 


The  second  term  on  the  right  is  just  u>  • I . The  first  term  on  the 
right  can  be  converted  to  an  expression  involving  H by  first 
writing  Hc  as  follows 


N 


H = I • u> 
c 


= -E 


R.  • w 

1C  1C 


i=l 


N 


V m.  R.  «v 

1 1C 

i=l 


where 


(A.  32) 


V 


• ou 


(A -33) 
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Hence 


Substituting  this  into  D now  yields 


D = 


7 <»< 


+ ou 


I + I 


I1  ) 


F rom  C 


J there  now  follows 


(A.  37) 


C 


H 

c 


J + I • J) 


(A.  38) 


Since  y in  equation  (A.  28)  is  block  diagonal  it  follows  that  in 
the  equation  G = y • G + K there  is  no  coupling  between  the 
rotational  and  translational  equations  of  motion.  Thus,  using  the 
y from  equation  (A.  28)  now  leads  to  the  following  equations  of  motion: 


H = C • H + L 
c c c 


(A. 39a) 


m = - c*  • UJ  + 7 • Lc  (A.  39b) 


Comparing  equations  (A.  39a)  and  (A.  la)  yields  the  requirement 


C • Hc  = 0 (A.  40a) 

j 


t 

1 

i 


I 
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Similarly,  comparing  equations  (A.  39b)  and  (A.  5a)  yields  the 
requirement 

I • C*  • uu  = UJ  • Hc  = H*  • uT  (A.  40b) 

Equations  (A.  40a)  can  easily  be  verified  with  the  use  of  C in 
equation  (A.  38): 

C • H = t ( ^ • H + H • J • H + I • a)*1.  J • H ) 

C C.  C C C c 

1 ~ ~ ~ -+ 

= y ( U)  • H + H • (D  + I . (IT  • ID  ) (A.  41) 

u c c 

= 0 


Note  that  the  first  term  on  the  right  of  this  equation  is  the  negative 
of  the  second  term;  the  third  term  is  zero  since  ao  • u;  = uj  x = 0. 
Thus  equation  (A.  40a)  is  verified.  Similarly 

C •U)  = ^-((JD*U)  + J»  . U)  + J . ou  . I . UJ  ) 

= J • ( H*  • uj+U).  Hc  ) (A.  42) 

— ~ t -* 

= J • H • w 
c 

Thus  equation  (A.  40b)  is  also  verified. 
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The  relationship  p = y.  a + a • yx  yields  the 
requirement 


C • I + I • C 


D + D 1 


(A. 43) 


But  from  equation  'A.  37)  it  follows  that 


D+Dt  = u>.I+I.tt)t  (A. 44) 

The  right  hand  side  of  equation  (A.  44)  is  known  to  be  equal  to  I . 
Hence  equation  (A.  43)  is  also  verified. 
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B.  INTERCONNECTING  TWO  RIGID  BODIES 


When  a system  of  bodies  is  interconnected  the  primitive 
system  equations  are  obtained  by  simply  stacking  all  the  individual 
bodies  equations.  Thus,  for  a two  body  system  the  equation 
G = n • g becomes 


M 


1 


(B.l) 


where  the  off-diagonal  terms  of  |i  are  zero  dyadics  O . Note 
that  superscipts  1 and  2 are  used  to  refer  to  bodies  1 and  2,  and 
subscripts  are  used  to  refer  to  the  center  of  mass  points  c^ 
and  c^.  The  momentum  formulation  equation  G - y • G = K 
is  given  by 


(B.2) 
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and  the  velocity  formulation  equation  p • g + V * • g = K is 
given  by 


— 

— 

— — 

— — 

— 

71 

& 

-i* 

c 

m1 

L1 

C1 

72 

ir2 

+ 

— >2t 

C 

L2 

c2 

~*i 

M 

V 

C1 

O 

V 

C1 

F 1 

— 2 
M 

V 

c2 

o 

— ► 

V 

C2 

F2 







__ 





(B.  3) 


B.  1.  Introduction  of  Relative  Velocities 

The  above  primitive  system  equations  have  12  degrees  of 

freedom.  But  in  a system  of  two  bodies  it  is  usually  convenient  to 

model  the  second  body  with  less  than  6 degrees  of  freedom.  In 

order  to  impose  restrictions  on  the  relative  motion  of  the  second 

body  it  is  convenient  to  first  transform  the  primitive  variables 

and  equations  to  new  variables  and  equations  so  that  the  impositions 

of  relative  motion  constraints  becomes  simpler. 

Let  h 2 be  a point  fixed  in  body  2 and  let  b ^ be  a point  fixed 

in  body  1.  The  position  vector  to  h_  from  b_  is  denoted  by  R,  , . 

1 L L n2  2 
Let  R.  , be  the  time  derivative  of  R,  , with  respect  to 

h2b2  h2b2  j. 

body  1.  Then  the  inertial  time  derivative  R^  ^ can  be  written  as 

2 2 
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X 


(B.  4a) 


R 


h2b2 


h2b2 


+ u> 


R 


h2b2 


In  order  to  put  bodies  1 and  2 on  a similar  footing,  let  be  a 
point  fixed  in  body  1 and  let  b^  be  at  the  inertial  reference  frame 
origin.  For  convenience,  this  origin  can  be  defined  by  an 
infinitesimal  body  0 with  an  angular  velocity  U)  ^ = 0 . Then 
analogous  to  equation  (B.  4a)  there  is  the  result 


Ri^  u.  = Rv  , + u>°  x R 

hlbl  hlbl  hlbl 


(B.  4b) 


—2  -»i 

Now  introduce  the  notation  V and  V as  follows: 


U = R 


h2b2 


(B.  5a) 


U = R 


hlbl 


= R 


hlbl 


= r,  = v. 


(B.  5b) 


The  additional  equalities  in  equation  (B.  5b)  follow  from  (B.  4b)  with 

A 


u>'°  = 0 


and  from  R^  ^ = r^  since  b^  is  at  the  inertial  origin. 


11  1X1 

The  symbol  r,  , of  course,  denotes  the  position  vector  to  the  point 
ni 


hj  from  the  inertial  origin;  the  inertial  time  derivative  of  r^  is 

v,  , the  linear  velocity  of  point  h.. 

1 1 

-*2  -*i 

The  relative  angular  velocities  O and  O are  now 
introduced  as  follows  : 


n2  = s;2  . i1 

n1  = ul1  . w°  = w1 


(B.  6a) 
(B. 6b) 
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The  velocity  transformation  equation  g = A • g now  takes 
the  form 


— — 

— 

— 

— — 

J1 

E 

d1 

w*2 

E 

E 

d2 

V 

C1 

R 1 h 
clhl 

E 

u1 

V 

C2 

Rt  . 

C2hl 

Rl  . 

C2h2 

E 

E 

d2 





— 

— — 

(B.  7) 


The  inverse  velocity  transformation  g = B • g is  given  by 


d1 

E 

S1 

-+? 

n 

-E  E 

u) 

= 

• 

— i 

u 

R , E 

V 

clhl 

ci 

-2 

U 

-R  . R . -E  E 

V 

clh2  C2h2 

c2 



— 

— 

(B.  8) 


The  elements  of  A and  B which  are  not  shown  are  zero  dyadics  O . 

The  equation  for  v and  v follow  from  taking  the  inertial  time 
C1  C2 

derivatives  of  the  equations 
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(B.  9a) 

(B.  9b) 


and  then  making  use  of  the  facts  that  R , and  R,  , are  fixed 

C1  1 b2  1 

in  body  1 whereas  R , is  fixed  in  body  2. 

c2  2 


Transformed  Variables  and  Equations 


Now  that  the  transformation  operator  A has  been  determined 
the  new  momentum  and  force  variables  can  be  determined  via  A^. 
Thus,  the  equation  G = A*  • G yields 


Now  introduce  the  following  notation 


H 

hl 

= H1 

C1 

+ 

R , 
clhl 

• P1 

»h 

hl 

= £2 
c2 

+ 

^C2hl 

• P2 

hl 

-5* 

hl 

+ 

»h 

hl 

H ^ 
h2 

= H2 
c2 

+ 

R , 
c2h2 

• P2 

P * 

= P 1 

+ 

P2 

(lj.Ha) 


(B! 


lib) 


(B.lllc) 


f 

(B.Jld) 


(B.  lie) 


For  i and  j equal  to  1 and  2,  is  the  angular  momentum  of 

body  i about  the  point  h.  . ^ is  the  angular  momentum 

1 -T 

of  the  system  of  two  bodies  about  the  point  h^,  and  P is  the 
linear  momentum  of  the  system  of  two  bodies.  Combining 
equations  (B.  10)  and  (B.  11)  now  yields 


G = 


*h 


1 


h; 


p2 


(B.  12) 
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\ 


The  equation  K = A • K yields  similarly 

r - 

r l 


K = 


(B.  13) 


where  is  the  torque  on  the  system  of  two  bodies  about  h^ 

-►  2 — 

is  ^he  torque  on  body  2 about  h^  , F is  the  force  on  the 

^ 1 ->2 

system  o$  two  bodies,  and  F is  the  force  on  body  2. 

* _ t 

Expanding  the  equation  M>  = A • p • A yields 


H = 


‘71 

1 hjhl 


h2hl 


1 

M ^crh 

,CI  1 


M2  b} 


^2hl 

» 


^lh2 


72 

h2h2 


M2  R4 

C2h2 


M2  R4  . 
C2h2 


M1  R 


M2  R 


cyhi 


C2h2 


M' 


M 


M 2 R 


C2hl 


M2  R 


C2h2 


2 

M 


M ' 


(B.  14) 


where  for  i,  j,  k equal  to  1 and  2, 
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*-**i  •*-►1  ■» 

^ = 1 + MR 


. • R U 

c.h.  c.h. 

1 J IK 


(B.  15a) 


T1  = T1 
hihi  Vi 


+ I 


hlhl 


(B.  15b) 


" l *— • l «-»  ? 

M = M + M 


M1  R = M1  R , + M2  R , 

CIhl  clhl  c2hl 


(B.  15c) 
(B.15d) 


I i u is  the  inertia  dyadic  of  the  system  of  two  bodies  about  the 
nlnl 

point  h^  , and  M is  the  mass  dyadic  for  the  system  of  two 
bodies.  The  point  cy  is  the  center  of  mass  of  the  system  of 

two  bodies,  and  R , is  the  position  vector  to  ct  from  h,  . 

c2  A 1 


v = 


Expanding  the  equation 

v = b • 

v • B*  yields 

J1 

-J1 

zj, 

hl 

-Z?, 

h2 

-J1 

U2  + 71) 

-zi 

hi 

+ zi  > 

2 n2 

hl 

W1 

hlhl 

-w1 

hlh2 

-1* 

(zjf  + Z £ ) 
h2  h2 

-w1 

h2hl 

<<2»2  + < 

2 2 


(B.  16) 
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where  for  i,  j,  k equal  to  1 and  2 


(B. 17a) 


+ ®c.h.  • ?i  • K.K  <B'17b| 

] k l j l Tc 

There  remains  the  evaluation  of  the  term  X = - y • G 

where  y = • y • B*  + A*  • B^  . When  the  reduced  form 

of  y is  used  (i.  e.  , when  X is  expressed  as  X = - yp.  • G ) 

G 

the  result  is 


(B.  18) 


0 


Consequently,  the  momentum  formulation  equation  of  motion, 
G + X = K,  takes  the  form 
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— — 
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~ -1 
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v,  • P 

h. 
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1 

1 

i 

^2 

~ ->2 

7 2 

Hh 

h2 

v,  • P 

h2 

h2 

+ 

= 

-7 

— “5 

p 

0 

F 

2 

0 

F 

— — 





(B.  19) 


This  result  is  well  known 

the  fundamental  result 

since  it  is  an  immediate 

consequence  of 

+ v X PA  = LA 

(B.  20a) 

a 

a 

a 

> 

n 

F A 

(B.  20b) 

where,  given  a system  A and  an  arbitrary  point  a,  H is  the 

—A  , a 

angular  momentum  of  A about  a,  L.  is  the  torque  on  A about  a, 

— a **■  —A 

P is  the  linear  momentum  of  A,  and  F is  the  force  on  A; 

v is  the  linear  velocity  of  the  point  a.  Equation  (B.  20)  is  used 
3.  _ 

twice  in  equation  (B.19):  once  with  A = 1 and  a = h^  , and  once 
with  A = 2 and  a = . In  the  first  case  system  A is  the  system 

of  two  bodies;  in  the  second  case  system  A is  body  2. 

Even  though  the  result  (B.19)  is  well  known,  equations  (B.14) 
for  p and  (B.  16)  for  v are  not  well  known.  In  fact,  these  equations 
have  only  been  published  in  connection  with  the  author's  works  on 
the  transformation  operator  approach  to  multi -body  dynamics 
[3],  [4],  [53]. 
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B.  3.  Imposition  of  Constraints 


Now  suppose  that  body  2 has  only  a rotational  degree  of  freedom 

— 2 -*2 
relative  to  body  1.  Then  U is  not  a free  variable;  i.  e.  , U 

-*2  - 

is  some  prescribed  function  of  time  (perhaps  U =0  and  the 

points  h?  and  b7  start  out  and  remain  coincident).  Therefore 
-2  ^ ^ — 2 

P is  not  a free  variable  either,  since  P can  be  determined 
algebraically  from  the  remaining  free  variables.  Consequently 
the  last  equation  of  equations  (B.  19)  can  be  removed  from  the  set. 

It  is  of  interest  to  relate  the  free  and  constrained  variables 

via  the  projectors  tt  1 = A.  • B.  for  i = 1 and  2.  The  free 

— 11  - —2 
velocity  is  and  the  constrained  velocity  is  ^ = U ; 


*1 


ft 


2 


(B.  21 ) 


Bj  is  now  determined  from 


g (B. 22a) 
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Similarly  is  determined  from 


(B.  22b) 


Now  g is  the  part  of  g due  to  g^  . Hence 


(B. 23a) 


A?  ' 


(B.  23b) 


From  equations  (B.23)  it  follows  that 


, where  v ' = v - U 

C2  C2 


(B.  24a) 


(B.  24b) 


Evidently,  g is  restricted  to  a 3-dimensional  subspace  of  the 
12-dimensional  space  which  contains  g. 

From  the  above  expression  for  and  B^  there  now 
follow  the  expression 


" = A1  • B1 


-R  , E 
c2h2 


(B.  25a) 
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o 


o 


o 


clh2 


R 


C2h2 


-E 


(B. 25b) 


Using  these  expressions  for  the  projectors,  it  is  readily  verified 
that 


g1  = n 1 • g (B. 26a) 

g2  = tt2.  g (B.  26b) 


Note  that  the  subspace  to  which  g1  is  restricted  is  evident  from  the 
form  of  the  projector  tt1  . 

The  free  force  is  K,  and  the  constraint  force  is  K-,  . For 

1 _ 1 2 L 
i = 1 and  2,  K is  the  part  of  K due  to  K.  . K and  K are  given 

as  follows 


„1  1 

K = TT 


K = 


R 


clh2 


R 


C2h2 


(B.  27a) 
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v2  2 
K = TT 


K = 


-R 


R 


clh2 


C2h2 


-F  ' 


F2 


F2 


(B. 28b) 


1 2 1 

From  these  expression  for  K and  K it  is  evident  that  K is 

2 

restricted  to  a 9-dimensional  subspace  whereas  K is  in  a 
3-dimensional  subspace. 
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